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Abstract 

An overview of the recent progress in the studies of nuclear mul- 
tifragmentation is presented. Special emphasis is put on the explo- 
ration of isotopic trends in nuclear multifragmentation and the possi- 
bilities to extract physical information related to the nuclear equation 
of state. Relevant experimental methods of isotope identification are 
described. The isotopic composition of fragments is used to extract 
the values of thermodynamical observables of the system undergo- 
ing multifragmentation such as temperature and chemical potentials. 
Various methods for extraction of thermodynamical variables are an- 
alyzed. An overview of methods of isotope thermometry, exploring 
the sensitivity of various yield ratios to temperature, is presented. 
An exponential scaling of relative isotopic yields from reactions with 
different neutron content, called isoscaling, is used to explore the evo- 
lution of the isospin degrees of freedom of the system. Finally, the 
nuclear equation of state and the isospin-asymmetric liquid-gas phase 
transition in the nuclear matter are discussed. 



Introduction 

Experimental studies, carried out during last two decades ( for review see 
and for update of the recent progress see e.g. [2] ), demonstrate that, in nu- 
clear reactions such as the nucleus-nucleus collisions in the Fermi energy do- 
main or the collisions of a high-energy light particle with a heavy target, the 
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amount of kinetic energy dissipated into the internal excitation energy is large 
enough to produce highly excited ( hot ) nuclei which may further de-excite 
via multifragmentation, a decay mode where the hot nucleus disintegrates 
into many fragments. In the nature, nuclear processes of such type can be 
observed in the interaction of cosmic rays with terrestrial and interplanetary 
matter. Such spallation reactions contribute to the observed isotopic compo- 
sition of the interstellar matter which can be used to study the properties of 
the cosmic rays and of their astrophysical sources ( for review see e.g. [3] ). 
With the advent of powerful accelerators of high-energy particles and, more 
recently, of heavy ions of intermediate and high energy, it became possible to 
study the process of multifragmentation in the laboratory. Extensive experi- 
mental studies of multifragmentation are conducted since eighties of the last 
century. As can be deduced from observed experimental data P3 |2] , during 
the process of multifragmentation the nuclear medium becomes hot, expands 
and is supposed to undergo a phase transition where two phases are created, 
a dense isospin-symmetric phase ( liquid ) and a dilute isospin-asymmetric 
phase ( vapor or gas ). Such a phase transition was predicted by theory ( for 
review see 0] ). The most common theoretical description of multifragmen- 
tation is based on the statistical model of multifragmentation [5], typically 
assuming a thermally equilibrated freeze-out configuration with hot ( excited 
) or cold fragments. The probabilities of various fragment partitions can be 
evaluated analytically using the grand-canonical approach. While the sta- 
tistical model of multifragmentation is quite successful in describing a wide 
range of fragment observables, it does not provide a specific information on 
the mechanism of the underlying process of fragment formation, apart from 
the assumption that at the freeze-out configuration a thermal equilibrium 
among the fragments is reached. The underlying process of cluster forma- 
tion can be related to the equation of state of asymmetric nuclear matter 
at sub-saturation densities, where the phase transition is supposed to occur. 
The coexistence line, typical for the one-component system such as symmet- 
ric nuclear matter, develops with introduction of isospin asymmetry into a 
complex coexistence region jK, typically suggesting a coexistence of the di- 
lute isospin-asymmetric phase and isospin-symmetric dense phase. Further 
theoretical investigations [7] suggest that instability modes responsible for 
the phase transition are typically a combinations of mechanical and chemical 
instability modes. 

In this article, we present an overview of the recent progress in the studies 
of nuclear multifragmentation. We place a special emphasis on the explo- 
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ration of isotopic trends in nuclear multifragmentation and the possibilities 
to extract physical information related to the nuclear equation of state. Ex- 
perimental studies of isotopic trends in nuclear multifragmentation have been 
made possible by use of the state-of-the-art detection techniques, which al- 
low to identify individual isotopes of light charged particles and intermediate 
mass fragments in wide angular ranges. The isotopic composition of frag- 
ments can be used to extract the values of thermodynamical observables of 
the system undergoing multifragmentation such as temperature and chemi- 
cal potentials ( free nucleon densities ). Various methods for extraction of 
temperature from isotopic composition are analyzed. An exponential scal- 
ing of the yields of specific isotopes from different reactions with the isotope 
neutron and proton numbers ( isoscaling ) is used to study the evolution of 
the isospin degrees of freedom of the system. Finally, the nuclear equation of 
state and the isospin-asymmetric liquid-gas phase transition in the nuclear 
matter are discussed. 




Multifragmentation theory and related topics 

Multifragmentation is a process where the final state consists of multiple 
fragments which are the remnants of a hot nucleus. The aim of the multi- 
fragmentation theory is to relate the properties of final fragments to prop- 
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erties of the original hot nucleus. Statistical phase-space models of mul- 
tifragmentation typically assume that the properties of final fragments are 
determined in the so-called "freeze-out" configuration, where the typical dis- 
tance of pre-fragments becomes, due to thermal expansion, larger than radius 
of the short-range nuclear interaction. It is, however, difficult to study the 
properties of the freeze-out configuration using the experimentally observed 
final fragments without additional model assumptions reflecting the reaction 
mechanism by which the hot nucleus was created. A detailed knowledge of 
the reaction mechanisms contributing to production of hot nuclei is neces- 
sary to properly characterize initial conditions of multifragmentation. The 
production of hot nuclei can be treated separately as a stage prior to multi- 
fragmentation, assuming that statistical equilibrium will be reached at later 
stage ( typically at freeze-out ). In a microscopic description of nuclear 
medium, multifragmentation can be related to a phase transition into dense 
and dilute phase. Isospin asymmetry of two phases may be different. Fur- 
thermore, multifragmentation can be treated dynamically using the transport 
theory. Such treatment offers the possibility to study the evolution of the 
nuclear medium during the whole process, including the conditions for the 
isospin-asymmetric phase transition. In this section we briefly review the 
most relevant aspects of the multifragmentation theory and related reaction 
theory. 

Statistical phase-space models of multifragmentation 

A main assumption of the statistical phase-space models describing various 
processes is that a certain stage of the process can be identified where the 
available phase-space determines the final probability distribution of exit 
channels. For instance, in statistical description of the fully equilibrated 
compound nucleus, the probabilities of various evaporation channels are de- 
termined by available phase space of the corresponding evaporation residue 
and the fission probability is determined by available phase space at the 
saddle-point on top of the fission barrier. In the case of simultaneous mul- 
tifragmentation, the phase-space is explored for the possible multi-fragment 
configurations ( partitions ), typically at the stage where short-range nuclear 
interaction freezes out and the identities of final fragments are determined. 
This is the so-called "freeze-out" configuration. The probabilities of various 
partitions can be determined statistically, exactly using the micro-canonical 
ensemble jS], or approximately using the canonical and grand-canonical en- 
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sembles jS]- Canonical and grand-canonical approximations typically reduce 
the mathematical complexity while raising questions on the applicability of 
such approximations since the nuclei are rather small systems. Detailed dis- 
cussion of the application of statistical mechanics to small systems can be 
found in 0. 

Here we briefly review the most commonly used phase-space model of 
multifragmentation, the Statistical Model of Multifragmentation ( SMM ) [5] . 
The model uses the grand-canonical approximation, the freeze-out volume 
depends on fragment multiplicity and internal excitation of fragments in 
the freeze-out configuration is considered, thus producing "hot" fragments. 
Analogous models with "cold" fragments also exist [S], typically the freeze- 
out volume is chosen approximately twice larger than in the case of hot 
fragments. Yield of fragments with neutron and proton numbers N and Z 
can be, within grand-canonical approach, characterized by formula 0: 



Y(N, Z) = Qnz^A 3 / 2 exp(-(l/T)(F NZ (T, V) - Nfi n - Zfi p )) (1) 
At 

where g^z is the ground state spin degeneracy, Vf is the free volume 
( Vf = V — Vq where V and V are the volumes at freeze-out and ground 
state ), Xt is the thermal wavelength, F^z is the intrinsic free energy of the 
fragment, fi n and /i p are the free neutron and proton chemical potentials, 
and T is the temperature. The free volume at freeze-out is determined as 
Vf = xVq- The parameter \ is determined as 

x = (l+ *(MVa -!))«- 1. (2 ) 

where M is the fragment multiplicity, d is nuclear diffuseness, A is the 
mass of the system and r is the nuclear radius parameter. It is essentially 
a product of the additional fragment surface area and nuclear diffuseness 
length. The values of x range from 0.2 to 2. The fragment free energy is 
at most part determined using the liquid drop model. The component of 
free energy corresponding to Coulomb interaction is determined using the 
Wigner-Seitz approximation. The Wigner-Seitz clusterization energy is cal- 
culated at volume V = (1 + n)V . The value of k can be related to average 
Coulomb barrier and thus is generally different from x, which controls the 
freeze-out of strong interaction. Typical value is k = 3.5 which is between 
the values for closely packed fragments ( k — 2 ) and for fission ( n — 5 ). The 
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statistical model of multifragmentation was used for extensive comparisons 
with experimental multifragmentation data, typically obtaining a consistent 
description of a wide range of experimental observables |o1 ITUl ITT] . 

As an alternative to models of simultaneous multifragmentation, the tra- 
ditional model of compound nucleus decay based on Hauser-Feshbach ap- 
proximation ^2] was extended to describe the emission of intermediate mass 
fragments ( IMF ) [JJJ] . The emission of such complex fragments is described 
as a binary split, essentially an asymmetric fission where the IMF is accom- 
panied by a heavy residue. The emission probability is determined by a 
value of the mass- asymmetric fission barrier height. The model of fission-like 
binary split describes reasonably well the onset of fragment emission, thus 
indicating that at the onset the complex fragment emission is a fission-like 
process jUj. Success of such approach prompted the attempts to describe 
multifragmentation using the model of Sequential Binary Decay ( SBD ). 
Several implementations ( e.g. the GEMINI code ^5] ) have been used ex- 
tensively for comparisons with experimental multifragmentation data, with 
reasonable success in description of e.g. inclusive mass and charge distribu- 
tions. However, the model of sequential binary decay fails to describe the 
results of the exclusive measurements at the zero angle JTj allowing to select 
a narrow region around the multifragmentation threshold while SMM per- 
forms rather well. An explanation is suggested by comparison of the model 
of sequential binary decay to transitional-state model of multifragmentation 
TIP where even for binary exit channel the number of available degrees of 
freedom is different. Furthermore, the model of sequential binary decay also 
fails to reproduce the observed fragment correlations |17j . 

Isospin-asymmetric liquid-gas phase transition 

Phase-space models of multifragmentation typically explore the properties 
of final fragments in the freeze-out configuration where nuclear interaction 
between the fragments does not play a role. The underlying mechanism 
leading to multifragmentation is nevertheless closely related to properties of 
the nuclear medium, which can be studied microscopically. Over the recent 
decades, thermal properties of the nuclear medium were studied theoretically 
using microscopic approaches such as temperature-dependent Hartree-Fock 
18J or relativistic mean-field [S| approximations. The microscopic theory of 
the nuclear matter predicts a phase transition of the liquid-gas type. Below 
a critical temperature, the two phases can coexist. The spatial distribution 
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of phases can be expected to develop from vapor bubbles in liquid at high 
density toward liquid droplets in vapor at low densities r J])j . A hot nucleus 
undergoing multifragmentation can be envisioned as rapidly expanding until 
it enters the spinodal region. At that stage multifragmentation occurs due 
to rapidly developing dynamical instabilities. When considering the isospin- 
asymmetric nuclear matter [UJ |2D] , the phase transition is supposed to lead 
to separation into a symmetric dense phase and asymmetric dilute phase. It 
has been discussed in literature (7^ to what extent such a phase transition is 
generated by fluctuations of density or concentration, typically suggesting a 
coupling of different instability modes. 
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Figure 2: A schematic description of the nuclear liquid-gas phase transition as a mech- 
anism of multifragmentation. 

The properties of clusters formed in the phase transition can be described 
by a variety of random cluster models. Simple model of clusterization was 
proposed by Fisher where the liquid drop expansion of the Gibbs free 
energy leads to the mass distribution of clusters 

P(A) oc A- T exp(-6(T)A 2/3 ) (3) 

where the surface term b(T) = 4irrQa(T)/T decreases monotonically to- 
ward zero as temperature approaches critical temperature T c . At critical 
temperature the mass distribution becomes a power law with critical expo- 
nent r. Similar properties are obtained using percolation models, lattice 
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models and cellular model of fragmentation ( an extensive overview of these 
models in the context of multifragmentation can be found in 0J, sections 
4,5. ). 

Calculations based on the microscopic nuclear theory [SJ 120] suggest that 
nuclear matter undergoes liquid-gas phase transition which is of the first or- 
der, the density difference between phases being the order parameter. Similar 
conclusions can be obtained also in the context of SMM (221 • However, the 
phase transition in the percolation models is of the second order and lat- 
tice models such as Ising model [2H] exhibit a complex phase diagram with 
both first and second order phase transitions. Furthermore, since theory of 
nuclear phase transitions is formulated for infinite systems while the experi- 
ments investigate multifragmentation of finite nuclei, it is necessary to assess 
the effect of finite size. As discussed in the critical behavior observed in 
random-cluster models can be related to scaling properties and self-similarity 
of the system at different scales, what suggests that observed properties of 
finite system represent the properties of infinite media. On the other hand, 
finiteness constraints of small systems can lead to a behavior analogous to 
criticality [3]. Obviously, for any experimental observation the effect of finite 
size should be assessed. 

Production of hot nuclei 

The statistical multifragmentation theory typically does not explicitly con- 
sider the evolution of the system prior to the freeze-out ( or underlying phase 
transition ). Usually it is difficult to determine the properties of the freeze- 
out configuration using the experimentally observed final fragments without 
additional model assumptions reflecting the reaction mechanism by which the 
hot nucleus was created. A detailed knowledge of the reaction mechanisms 
contributing to production of hot nuclei is necessary to properly characterize 
initial conditions of multifragmentation and to verify the model assumptions 
used by multifragmentation theory. 

In order to undergo multifragmentation, intrinsic excitation energy of the 
hot massive nuclei should exceed 2-3 AMeV. Such excitation energies can 
be reached in nucleus-nucleus collisions with the nuclear beams of interme- 
diate energies ( 20 - 100 AMeV ). In collisions of a high-energy light particle 
with a heavy target the necessary energies range from hundreds of AMeV to 
AGeV. These are the reaction domains which have been used extensively for 
multifragmentation studies. 



8 



1 



Figure 3: Development of the reaction mechanisms in nucleus-nucleus collisions in the 
Fermi energy domain. Cross section contributions are plotted schematically as a function 
of initial angular momentum. 

A large variety of reaction mechanisms has been observed in nucleus- 
nucleus collisions in the Fermi energy domain ( 20 - 50 AMeV ) depending 
on the impact parameter, projectile-target asymmetry and the projectile en- 
ergy. The reaction mechanisms typically observed are |24j : 

- peripheral elastic and quasi-elastic ( QE ) scattering/transfer reactions 
around the grazing impact parameter. 

- deep inelastic transfer ( DIT ) reactions at semi-peripheral impact pa- 
rameters with partial overlap of the projectile and target and a significant 
part of the relative kinetic energy transferred into internal excitation energies 
of the projectile and target. 

- incomplete fusion ( ICF ) reactions at central impact parameters with a 
typical participant-spectator scenario. At energies around the Fermi energy 
the participant zone typically fuses with one of the spectator zones (typically 
heavier) thus creating a highly excited composite nucleus. At most central 
impact parameters the complete fusion ( CF ) can occur. 

- pre-equilibrium emission of direct particles, caused by onset of two-body 
nucleon-nucleon collisions at central and mid-peripheral impact parameters. 
Pre-equilibrium emission typically precedes the ICF/CF and DIT reactions. 
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While such a scheme is rather well justified by experimental data 
there are still many open questions related to the reaction mechanism, specif- 
ically concerning the development of isospin-asymmetry in nucleus-nucleus 
collisions in the Fermi energy domain. Studies of isospin degrees of freedom 
in the nucleus-nucleus collisions in the Fermi energy domain reveal many 
interesting details of the reaction scenario. Recently, an enhancement in the 
production of neutron-rich nuclei was observed in peripheral nucleus-nucleus 
collisions [22], which can be related to the effect of neutron-rich surface of 
the target. 

At projectile energies of hundreds of AMeV and above the impinging light 
particle initiates an intra-nuclear cascade [2H] of nucleon-nucleon collisions 
which is the main mechanism leading to production of hot nuclei. The heavy 
remnant after the intra-nuclear cascade is highly excited and typically un- 
dergoes multifragmentation at wide range of impact parameters fSj- In the 
nucleus-nucleus collisions at similar energies the participant zone develops in- 
dependently of two spectator zones and typically disintegrates due to intense 
cascade of nucleon-nucleon collisions, with possibility to produce sub-nuclear 
particles via excitation of hadronic degrees of freedom. The spectator zones 
are relatively colder, nevertheless typically hot enough to undergo multifrag- 
mentation [27]. 

Dynamical description of multifragmentation 

As an alternative to two-step reaction scenario, where reaction dynamics 
is attributed to the early stage of nucleus-nucleus collision while multifrag- 
mentation of the hot nucleus is described statistically, the whole process 
can be described in unified fashion by dynamical theory with explicit time- 
dependence. 

In the microscopic time-dependent Hartree-Fock theory [2H], an exten- 
sion of the Hartree-Fock theory of the ground state, the nucleus-nucleus 
collision is represented by time evolution of the one-body mean field, which 
defines the collective degrees of freedom of the system. The time-dependent 
Hartree-Fock theory does exhibit collective behavior which reminds many 
experimentally observed phenomena |2B], it nevertheless does not include 
two-body dissipation via collisions of nucleons which is supposed to play 
significant role at energy regimes where the hot nuclei are produced. 

Two-body interaction can be included in the framework of quantum ki- 
netic theory as represented by the Boltzmann transport equation, which de- 
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Figure 4: Angular coverage of the multidetector array FAUST 

scribes the reaction dynamics in terms of both mean-field and two-body 
dissipation. Various semi-classical methods featuring the concept of pseudo- 
particle are used to simulate the solution of the Boltzmann transport equa- 
tion. The most common approximations are Vlasov-Uhling-Uhlenbeck ( 
VUU ) [H], Boltzmann-Uhling-Uhlenbeck ( BUU ) Landau- Vlasov ( LV 
) [HI], Boltzmann- Nordheim- Vlasov ( BNV ) jH2] and Quantum Molecular 
Dynamics ( QMD ) methods. Results of such simulations typically offer 
a qualitative description of the main features of reaction dynamics including 
multifragmentation. Using the transport theory, the properties of the nuclear 
equation of state such as nuclear incompressibility and symmetry energy can 
be investigated phenomenologically. 

Experimental aspects of multifragmentation 
studies 

Experimental studies of isotopic trends in nuclear multifragmentation were 
enabled by the recent progress in acceleration and detection techniques. 
A crucial condition is unambiguous isotopic identification of detected frag- 
ments, which can be achieved using high-resolution particle telescopes. In 
order to reconstruct the multifragmenting system and to understand the 
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Figure 5: Diagram of angular coverage of detector rings in the Att array NIMROD |35) . 

production mechanism of the hot nuclei it is desirable to detect the charged 
particles in a wide angular range, ideally in the 47r geometry. In this section 
we briefly review the experimental devices ( using the facilities installed at 
Texas A&M University as examples ) and essential methods relevant to the 
studies of isospin degrees of freedom in nuclear multifragmentation. 

Arrays of charged particle detectors with large angular 
coverage 

During the process of multifragmentation, multiple charged particles are 
emitted in a wide angular range. Angular distribution of fragments in the 
laboratory frame is determined by the production mechanism of the hot nu- 
cleus and by the choice of the projectile and target nuclei. The design of 
arrays of charged particle detectors must be optimized in order to detect 
and resolve emitted charged particles from the expected multifragmentation 
source. This can be achieved by choice of angular coverage, granularity and 
detection thresholds. 

A state-of-the-art forward array, the multidetector array FAUST |34| . in- 
stalled at the Cyclotron Institute of Texas A&M University, covers forward 
angles between 2 - 35 ° where the fragments originating from the projectile- 
like source can be expected. It consists of 68 charged particle telescopes 
arranged into 5 rings, as shown in Fig. 0] Each particle telescope consists of 
300 fim thick silicon detector followed by 3 cm thick CsI(Tl) crystal. Angular 
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Figure 6: An example of the raw Si-Si spectra The 4 He , 9 Be and 15 N lines used 
for calibration are shown as red lines. The data from reaction 124 Sn+ 124 Sn at 28 AMeV 
were collected in ring 8 of NIMROD. 

ranges covered by the rings are chosen in order to distribute the multiplic- 
ity of detected particles evenly, thus avoiding the cases when one telescope 
is hit by multiple particles. The mass and atomic number of the detected 
charged particles can be identified up to Z=6. Forward arrays such as FAUST 
are optimized for the study of projectile multifragmentation, specifically for 
asymmetric reactions where target nucleus is much heavier than the projec- 
tile. 

In order to study symmetric reactions where a hot composite source is 
created it is necessary to cover largest possible angular range, ideally Air. 
For instance, the multidetector NIMROD [33] at the Cyclotron Institute of 
Texas A&M University provides such a coverage using 96 detector modules 
arranged in 12 rings ( see Fig. |3] ). Light charged particles ( Z<2 ) are 
detected and identified in the whole range while a subset of modules ( 2 per 
ring ) offer isotopic identification up to Z=8. Furthermore, atomic number 
of heavy residues can be resolved using ionization chambers and neutron 
multiplicity can also be determined since the array is surrounded by a Air 
neutron multiplicity counter ( neutron ball.). 
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Figure 7: Calibrated Si-Si spectrum from Fig. The isotopic lines used for identification 
are shown. 

Identification of isotopes 

An important step in the off-line analysis of data from the multidetector ar- 
rays with large angular coverage such as FAUST and NIMROD is to identify 
as many fragment species as possible. The method of isotope identification 
is based on well known particle telescope technique in which the isotopes are 
resolved in two-dimensional AE-E spectra of energy losses in two detectors 
( typically a thin one followed by a thick one ). We developed and success- 
fully used [HI)] a method which enables to carry out identification and energy 
calibrations simultaneously using a minimization procedure applied to two- 
dimensional spectra. In the experimental spectra, the lines for three known 
isotopes ( typically the most characteristic isotopes such as 1 H, 4 He, 9 Be ) are 
assigned and calibration is carried out by a minimization procedure where 
these lines are fitted to corresponding calculated energy losses for a given 
AE-E telescope. The calibrations coefficients are thus obtained as optimum 
values of minimization parameters. 

For Si-Si calibration, the formula E = a + a\X + a 2 \/X~ was used for 
both Si detectors ( X being the amplitude of electronic signal ). Additional 
constraints were applied in the minimization functional in order to position 
correctly the punch-through points. Figure 10 shows an example of the ex- 
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Figure 8: An example of the raw Si-Csl spectra with de-calibrated isotopic dependences 
of energy losses. The 4 He , 7 Li and 9 Be lines were used in calibration. The data from 
reaction 40 Ca+ 27 Al at 45 AMeV were collected in ring E of FAUST. 

perimental spectra with selected 4 He , 9 Be and 15 N lines and the calibrated 
spectrum with calculated AE-E lines for a wide range of isotopes is shown in 
Fig. [7| The data used were collected in ring 8 of the NIMROD array in the 
reaction 124 Sn+ 124 Sn at 28 AMeV. One can see that the overall agreement in 
the shape of calculated and calibrated energy losses is very good, both inter- 
polating within and extrapolating beyond the range of selected isotopes ( at 
least to the neighboring element ). The way to identification of different iso- 
topes by linearization of experimental spectra and comparison to calculated 
lines is straightforward. 

For Si-CsI(Tl) calibration a similar procedure can be carried out. The 
difficulty when compared to Si-Si telescope is caused by the fact that the 
measured signal in CsI(Tl) crystal is the light output which depends on 
particle energy non-linearly ( including explicit dependence on particle mass 
and charge ). An approximate formula of Tassan-Got [37] for the energy 
calibration of CsI(Tl) detector was used 
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Figure 9: Isotope resolution obtained after linearization of Si-Si spectra. Data from the 
reaction 124 Sn+ 124 Sn at 28 AMeV were collected in the ring 8 of NIMROD [551155] . 

where L is the light output and p = rjAZ 2 where rj is a calibration parameter. 
An example of obtained agreement of experimental and calculated lines is 
given in Figure [HI using the data detected in the FAUST array in the reaction 
40 Ca+ 27 Al at 45 AMeV. The calibration of Si detector was obtained using 
radioactive source of a-particles and the calibration coefficients for CsI(Tl) 
detector were obtained using minimization procedure. Due to explicit mass 
and charge dependence the calibration must be preceded by identification, 
so in this case the calculated AE-E lines are de-calibrated using calibration 
coefficients and projected onto experimental spectrum. Again the overall 
agreement is quite good for both interpolation and extrapolation. 

The examples of obtained isotope resolution ( after linearization ) are 
given in Figs. 191101 Data from the reaction 124 Sn+ 124 Sn at 28 AMeV were 
collected in the ring 8 of NIMROD. One isotope line for each Z was used in 
linearization. The examples of particle spectra calibrated using the above 
methods are given in Fig. ^2 Data used in Fig. ^2 were collected in the 
reaction 124 Sn+ 28 Si at 28 AMeV using NIMROD j^HlEi- The experimental 
spectra are compared to simulations [21] described in previous section. 
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Figure 10: Isotope resolution obtained after linearization of Si-Csl spectra. Data from 
the reaction 124 Sn+ 124 Sn at 28 AMeV were collected in the ring 8 of NIMROD [5511591 . 

Event characterization 

Modern multidetector arrays with large angular coverage offer a possibility 
to study not only the inclusive observables such as isotope yields and parti- 
cle spectra but allow to study development of the fragment properties with 
evolving reaction dynamics. Typical goal of the exclusive study is to iden- 
tify dynamical properties of the hot multifragmentation source being created 
in the early stage of the nucleus-nucleus collision. Here we provide some 
examples for both peripheral and violent collisions. 

Using the forward arrays like FAUST one typically observes the mul- 
tifragmentation of the hot quasiprojectile which originates from peripheral 
and mid-peripheral nucleus-nucleus collision. We present here an example 
of the projectile multifragmentation in the reaction of 28 Si beam with 112 Sn 
and 124 Sn targets [10] - A hot multifragmenting source was identified us- 
ing the calorimetry technique for the events where all emitted fragments 
were isotopically identified. The velocity distributions of the reconstructed 
quasiprojectiles with total charge Z to t = 12 — 15 were determined. Resulting 
velocity distributions for the projectile energies 30 and 50 AMeV are given 
in Fig. iriZkb. Solid squares represent the reactions with 112 Sn target and 
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open squares reactions with 124 Sn. For a given projectile energy, the mean 
velocities and widths of distributions are practically identical for both tar- 
gets. The shapes of velocity distributions are close to Gaussian. Velocity 
distributions are symmetric and no significant low or high energy tails are 
observed. Thus, the reconstructed quasiprojectiles appear to originate from 
a single source which can be identified with the quasiprojectile. 

An apparent charged particle excitation energy of the quasiprojectile can 
be reconstructed for each projectile fragmentation event from the energy 
balance in the center of mass frame of the quasiprojectile 

K PP = Y,( T ? P + Am f)- Am QP (5) 
/ 

where T® P is the kinetic energy of the fragment in the quasiprojectile 
frame and Am/ and Ajtiqp are mass excesses of the fragment and quasipro- 
jectile, respectively. Even if emitted neutrons are not included in this ob- 
servable, it is useful for the light fragmenting systems, where the neutron 
emission is not expected to dominate. The distributions of the apparent 
quasiprojectile excitation energies reconstructed from fully isotopically re- 
solved events are shown in Fig. US] for both Z tot = 14 and Z tot = 12 — 15 ( 
solid circles and squares respectively ) along with the distributions obtained 
using DIT/SMM simulation ( solid histograms ). The agreement of the 
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Figure 12: Experimental velocity distributions of the fully isotopically re- 
solved quasiprojectiles with Z tot = 12 — 15 [10 ( solid and open squares mean 
112 Sn and 124 Sn target, (a) and (b) the projectile energy 30 and 50 AMeV, 
respectively ). Solid lines mean Gaussian fits. 



simulated and experimental apparent quasiprojectile excitation energy dis- 
tributions is quite good. The onset of multifragmentation into channels with 
Zf < 5 in the low energy part is described with good precision for all sets of 
data. 

While it is possible to fully reconstruct the projectile-like sources in the 
mass range A=20-30 in peripheral collisions, with increasing violence of the 
nucleus-nucleus collision a hot nucleus is formed from the parts of both the 
projectile and target. The excitation energy of such a source is high enough 
for emission of several intermediate mass fragments ( IMF ). Information 
on both mass and charge of emitted IMFs can be crucial to determine the 
emission mechanism and thus properties of the emitting source. Various 
techniques have been developed in order to characterize the impact parameter 
of the collision ( for a detailed study see [40J ). A simple method which is 
nevertheless widely used to determine the centrality of collision is based on 
the assumption that the particle multiplicity increases monotonously with 
impact parameter [IT] . However, such assumption may be too simplistic and 
more sophisticated methods using the combinations of several observables 
are developed jlUj- As an example, we present here a method [3Bj where 
the centrality bins are selected in the two-dimensional histogram of IMF 
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Figure 13: Distributions of the reconstructed apparent excitation energies of 
the quasiprojectiles ^Uj. Symbols mean experimental distributions of the set 
of fully isotopically resolved quasiprojectiles with Z tot = 14 ( solid circles ) 
and Z tot = 12 — 15 ( solid squares ). Solid histograms mean simulated distri- 
butions for corresponding data sets, (a) - 28 Si(30 AMeV)+ 112 Sn, (b) - 28 Si(30 
AMeV)+ 124 Sn, (c) - 28 Si(50 AMeV)+ 112 Sn, (d) - 28 Si(50 AMeV)+ 124 Sn. 



multiplicities ( Mimf ) vs multiplicities of charged particles ( Mcp ) by 
parallel lines chosen so that the most central bin integrates the events with 
highest multiplicities of charged particles and highest IMF multiplicities ( see 
Fig. ITU ) ) with the intermediate event classes, thus reflecting the effect of 
finite size. Fig. IT4*k shows that the shapes of spectra develop with increasing 
centrality. 

Event shape analysis was carried out for two-dimensional Mjmf vs. Mcp 
cuts as well as for the one-dimensional M IMF and Mcp cuts. All types of 
centrality cuts exhibit similar trends in terms of sphericity and coplanarity ( 
see Fig. nevertheless some differences can be observed in the most cen- 
tral bins. When considering the total number of events selected by various 
criteria, the statistics in the most central bin of the two-dimensional selection 
increased by two orders of magnitude compared to the traditional M C p se- 
lection ( see Fig. H5l). Such a gain in statistics is caused by better separation 
of central and peripheral events. The one- dimensional cuts typically admix 
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Figure 14: (a) - Inclusive spectrum of a-particles ( top ) from reaction 
124 Xe+ 124 Sn and spectra for 6 centrality bins ordered from top to bottom 
by increasing centrality [SB]- (b) - Centrality cuts in the M IMF vs. M C p 
histogram. Skewed lines indicate centrality cuts used in (a), horizontal and 
vertical lines indicate one-dimensional cuts in Mjmf and Mcp histograms. 



high-multiplicity peripheral events to a smaller number of genuine central 
events, thus decreasing the overall centrality. The use of two-dimensional 
criteria allows more detailed studies of central events. 



Isotopic trends 

The isotopic composition of the multifragmentation products such as light 
charged particles ( LCP ), intermediate mass fragments (IMF ) and heavy 
residues is an experimental information which can be used to trace the dy- 
namics of isospin degrees of freedom during the multifragmentation and, 
using grand-canonical picture, to extract the values of thermodynamical ob- 
servables of the system undergoing multifragmentation such as temperature 
and free nucleon densities. The isotopic composition can be studied both as 
a bulk or in terms of yield ratios of specific isotopes. A variation of the neu- 
tron content of reaction products can be used to characterize the response of 
the isotopic composition to variation of initial conditions. 
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Figure 15: Test of centrality selections in terms of sphericity (a), coplanarity 
(b), event statistics (c) and flow angle (d) for two-dimensional M IMF vs. Mqp 
cuts ( triangles ) as well as for the one- dimensional M IMF ( squares ) and 
Mqp ( circles ) cuts 

Isotopic distributions 

A wide variety of fragment species is typically produced in multifragmenta- 
tion of heavy nuclei. In order to obtain information on isospin dynamics, an 
obvious choice is to study the behavior of isotopic chains. As an example we 
provide in Fig. ^] isotopic distributions of fragments with Z=2-8 from four 
reactions of 124 Sn, 124 Xe beams with 112 ' 124 Sn targets at 28 A MeV [Ml E3 
studied at the Cyclotron Institute of Texas A&M University using the An 
multi-detector array NIMROD ..35j. The isotopic distributions from different 
reactions are distinguished by different symbols. Each panel shows the iso- 
topic distributions measured at three angular ranges corresponding to rings 
2, 5 and 8 ( from top to bottom ). 

As one can see, the overall shapes remain quite similar at different angular 
ranges. It is interesting to note that the isotopes with maximum yield for 
Z>3 are typically those with A=2Z+1 for both even and odd Z. The effect 
of different reactions is most pronounced at the proton and neutron-rich 
tails of distributions, neutron-rich isotopes being more abundant in neutron- 
rich systems and vice versa, while the central part corresponding to /3-stable 
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Figure 16: Isotopic distributions of fragments with Z=2-8 from four reactions 
of 124 Sn, 124 Xe beams with 112 < 124 Sn targets at 28 AMeV PUSH!, measured 
using the 4tt multi-detector array NIMROD [33] at three angular ranges. 



nuclei is relatively insensitive. In any case, such comparisons of the shapes of 
isotopic distributions are rather qualitative. A quantitative analysis can be 
carried out in terms of integral observables such as momenta of distributions. 

To express the effect of experimental conditions such as projectile-target 
combination or detection angle quantitatively one can evaluate the mean 
N/Z- ( neutron to proton ) ratio for each isotopic chain. In Fig. El are 
given the mean N/Z-ratios of emitted LCPs and IMFs with Z=2-8 from four 
reactions obtained at different angles 38j. The overall dependence is simi- 
lar at all angular ranges. He-isotopes are dominated by a-particles. Mean 
N/Z ratios are highest for Li- and Be-isotopes and decrease gradually with 
increasing atomic number of IMFs. When comparing mean N/Z ratios from 
four reactions two characteristic patterns can be distinguished. At most for- 
ward angles the mean N/Z ratios of Li- and Be-isotopes appear to track with 
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Figure 17: Mean N/Z-ratios of emitted LCPs and IMFs with Z=2-8 from 
four reactions obtained at different angles [SHI EH] ■ 



the isospin asymmetry of the composite projectile-target system what leads 
to grouping of experimental point into a typical 1-2-1 pattern. On the other 
hand the isotopes with Z=5-8 appear to track with the isospin asymmetry of 
the target nucleus what creates a 2-2 pattern. At angles 11.0-14.9° the 1-2-1 
pattern can be identified for Li- and Be-isotopes while for heavier fragments 
the patterns vary. At angles 28.6-35.5° one can recognize the 2-2 pattern 
for Li-isotopes while at heavier fragments the 1-2-1 pattern dominates. In 
principle, one can identify the 1-2-1 pattern with the fragments emitted from 
a hot composite source and 2-2 pattern to the emission from a quasiprojec- 
tile. The fragments with Z=5-8 at forward angles possibly originate from 
the binary de-excitation channels of the quasiprojectile while the 2-2 pattern 
of Li-isotopes at the most central angles suggests either backward emission 
from the quasiprojectile or emission from quasitarget. Similar isospin asym- 
metry patterns can be identified also in the isotopic yield distributions in 
Fig. Thus, the mean values of N/Z ratio appear to exhibit an inter- 
esting information on isospin dynamics. The relative variation of the mean 
N/Z values is nevertheless rather small when compared to the variation of 
isospin-asymmetry of the projectile-target combination and more sensitive 
isospin tracers should be identified. 
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Yield ratios 

As demonstrated in the previous subsection, variation of the projectile and 
target nuclei influences the integral characteristics of isotopic distributions 
such as mean values of N/Z ratio rather weakly. Nevertheless, when observ- 
ing the evolution of isotopic distributions in different reactions and angular 
ranges, one can notice that the relative yields of neighboring isotopes are 
rather sensitive ( especially at the tails of distributions ). Such yield ratios 
have been widely used in experimental studies of multifragmentation data 
during last 15 years 
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Figure 18: Three isotopic yield ratios for the multifragmentation data from 
four reactions of 124 Sn, 124 Xe beams with 112 ' 124 Sn targets at 28 AMeV [HBJ 
at three different angular ranges. The data from different reactions are 
characterized by a N/Z-ratio of the composite system ( N/Z c s . ). 

As an example, we present in Fig. UHlthe values of three isotopic yield 
ratios, again for the multifragmentation data jHHl EH! from four reactions of 
124 Sn, 124 Xe beams with 112 ' 124 Sn at 28 AMeV. In most cases one can observe 
an approximately exponential dependence of the isotopic yield ratio on the 
N/Z ratio of the composite system ( N/Z c s . ), which can be expected accord- 
ing to formula when assuming linear dependence of neutron chemical 
potential on N/Z c . s .. The yield ratios depend rather strongly on the angular 
range. Such a dependence can be caused by different isospin dynamics at 
different angles. 
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Figure 19: Three isobaric yield ratios for the multifragmentation data from 
four reactions of 124 Sn, 124 Xe beams with 112 - 124 Sn targets at 28 AMeV [HBJ 
I3TI] at three different angular ranges. The data from different reactions are 
characterized by a N/Z-ratio of the composite system ( N/Z c s . ). 



Furthermore, we present in Fig. flQlthe values of three isobaric yield ra- 
tios for the same multifragmentation data. Specifically, the isobaric ratios 
are evaluated for pairs of mirror nuclei. The sensitivity of isobaric ratios to 
N/Z c s is significantly stronger than in the previous case, especially for the 
ratios 9 Be/ 9 B and U B/ 11 C where the fragments 9 B and 11 C are populated 
rather weakly. The overall behavior is nevertheless rather irregular. Espe- 
cially the observed pattern at most forward angles suggests a contribution 
from the cold process such as projectile break-up, which appears to con- 
tribute considerably to production of isospin-asymmetric species at forward 
angles. Such observation is in agreement with available experimental data 

A behavior of the yield ratio of mirror nuclei 3 H/ 3 He was investigated in 
the quasiprojectile multifragmentation of a 28 Si beam with 112 Sn and 124 Sn 
targets at 30 and 50 AMeV j33]. As described in experimental section, the 
properties of the system which actually undergoes multifragmentation, in- 
cluding N/Z ratio ( N/Zqp ), are known with good precision. Fig. |2*U1 shows 
dependence of the isobaric yield ratio Y( 3 H)/Y( 3 He) on N/Zqp for nine bins 
of the apparent excitation energy per mass unit of the quasiprojectile ( e* app 
). The experimental data are represented as squares and the lines represent 
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Figure 20: Dependence of the yield ratio Y( 3 H) /Y( 3 He) on N/Z ratio of the 
isotopically resolved quasiprojectiles with Zq P = 12 — 15 for nine bins of e* pp 

jSj. 



the exponential fits. The observed exponential dependence in all excitation 
energy bins is in agreement with formula (JTJ). The slopes are steepest at low 
excitation energies and become flatter with increasing excitation energy. 

The apparent compliance of the yield ratio behavior to the grand-canonical 
picture as represented by the formula (JTJ) suggests the possibility to extract 
the values of thermodynamic observables, in particular temperature. The 
procedure for extraction of the temperature was suggested by Albergo 
When expressing the isotopic yield within the equilibrium limit of the grand- 
canonical ensemble, a simplified analogue of the formula (|TJ is obtained 



Y(N, Z) = F(V, T)g(N, Z) exp{(l/T)(B(N, Z) + N^ n + Z^)) (6) 

The isotopic yield is thus essentially governed by free neutron and pro- 
ton chemical potentials fi n , p , temperature T and fragment binding energy 
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B(N,Z). The factor F(V,T) encompasses information on the global prop- 
erties of the fragment partition. The ground state spin degeneracy factor 
g(N, Z) represents the internal state sum of the fragment, which is restricted 
to ground state only. In oder to extract temperature, a double isotope ratio 
is constructed [IS] 



Y(Ni + 1, Z 1 )/Y(N 1 , Zt) g{N x + 1, Z 1 )/g(N 1 , Z t ) 



x exp 



Y(N 2 + 1, Z 2 )/Y(N 2 , Z 2 ) g(N 2 + 1, Z 2 )/g(N 2 , Z 2) 
.BjNi + 1, gi) - B(N U gi) - B(N 2 + 1, Z 2 ) + B(N 2 , Z 2 

T 



x 



and the temperature T can be obtained as 



AB 



1234 



, Y(N 1 +1,Z 1 )/Y(N 1 ,Z 1 ) 
^ u Y(N 2 +l,Z2)/Y(N 2 ,Z2) 



■) (7) 



(8) 



where a = §|tSt^tj and = B ^ + 1. ^ - B ^ ^ ' 

B(N 2 + 1, Z 2 ) + B(N 2 , Z 2 ). Thus an estimate of the grand-canonical temper- 
ature can be determined using the isotope yields and some of their ground 
state characteristics. 
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Figure 21: Systematics of corrected isotopic ratios Y N>Z /Y N+ i^ z from four 
reactions plotted against the difference of binding energies 



The formula (JEJ) was used extensively in experimental studies during the 
last decade |2Z1 HZ|- Typically, the temperature determined using various 
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thermometers ( sets of isotope yields ) exhibited large spread, possibly due 
to side-feeding from secondary decay. A method of correction for secondary 
decay was proposed by Tsang et al. 48]. In [3Bj we employed a graphical 
method to extract the average temperature. The isotopic yield ratios, cor- 
rected for ground state spin degeneracy, are plotted as a function of binding 
energy difference ( see Fig. |2H)- The data from four reactions of 124 Sn, 124 Xe 
beams with 112,124 Sn targets at 28 A MeV, used already in Figs. 1181191 are in- 
vestigated. Individual double isotope ratio thermometers can be represented 
( in the logarithmic scale ) when connecting the two points corresponding to 
used isotope ratios by a straight line. The inverse of the slope of such line 
is then the temperature. When assuming that the effect of side-feeding in 
the secondary decay is distributed randomly, an average temperature can be 
obtained by a linear fit of the systematics of all isotopic yield ratios. 

An information on the temperature can be obtained also from the slope 
of dependence of the isobaric yield ratio Y( 3 H) /Y( 3 He) on neutron to proton 
ratio of the quasiprojectile N/Zgp. The data from reactions of a 28 Si beam 
with 112 Sn and 124 Sn targets at 30 and 50 A MeV ( see Fig. [2"U1 ) were used 
|44| and the extracted temperatures compared well with the temperatures 
obtained using the double isotope ratio thermometer ( 2 H/ 3 H)/( 3 He/ 4 He). 
An assumption was employed that the values of chemical potentials can be 
approximated by the values of separation energies. The reason why such 
approximation is applicable for hot expanded nucleus can be in the possible 
interplay of the effect of expansion ( reducing absolute values of chemical 
potentials ) and separation of the isospin-asymmetric free nucleon gas from 
the isosp in-symmetric heavier liquid phase [44], as indicated by observed 
inhomogeneous isospin distribution among light charged particles and heavier 
fragments |4*9*] . 

The isotopic yield can be, within the equilibrium limit of the grand- 
canonical ensemble, expressed in terms of the density of the free nucleon gas 



where p ntP oc exp(/i n p /T). The values of yield ratios of mirror nuclei with 
N — Z = 1 then can be expressed as 
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Y(N, Z) ex Vp»p*g(N, Z) exp(S(iV, Z)/T) 



(9) 



Y(N,Z) g(N,Z) p n 
Y(Z,N) g(Z,N)p p 



exp(AB/T) 



(10) 
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thus relating the yield ratio to the density ratio of the free neutron and 
proton gas. The yield ratios of the mirror nuclei 3 H/ 3 He, 7 Li/ 7 Be and 
u B/ n C were used in work [50J to determine the ratio p n / p p by an exponential 
extrapolation toward AB = 0, yielding the values approximately 2-5. While 
such approximation is valid for grand-canonical ensemble of infinite systems, 
when dealing with the finite systems such as nuclei the values of free nucleon 
chemical potentials depend on the fragment partition in a way guaranteeing 
conservation of the mass and charge of the system. Thus the value of the 
ratio p n / p P may differ strongly for different event classes and no globally 
applicable value may exist in the finite system. Furthermore, the effect of 
secondary emission may also contribute to uncertainty in determination of 

Pn/Pp- 

Isoscaling 

As shown in the previous subsection, the fragment yield ratios can be used 
to extract thermodynamical observables of the fragmenting system such as 
temperature and chemical potential. In the context of isotopic distributions, 
the fragment yield ratios represent the details of the distribution sensitive 
to isospin degrees of freedom. Similar sensitivity can be explored globally 
by investigating the ratio of isotopic yields from two processes with differ- 
ent isospin asymmetry, essentially dividing the two isotopic distributions in 
point-by-point fashion. When employing the formula such a ratio will 
depend on N and Z as follows jSI] 

R 21 (N,Z) = Y 2 (N,Z)/Y 1 (N,Z) = Cexp(aN + ftZ) (11) 

where a = Ap n /T and ft = Ap p /T, Ap n and Ap p are the differences in 
the free neutron and proton chemical potentials of the fragmenting systems. 
C is an overall normalization constant. Alternatively the dependence 
can be expressed as 

R 21 (N, Z) = Y 2 {N, Z)/Yl(N, Z) = Cexp(a'A + ft'{N - Z)) (12) 

thus introducing the parameters which can be related to the isoscalar and 
isovector components of free nucleon chemical potential since a' = A(p n + 
/x p )/2T and ft' = A(p n - p p )/2T. 

An exponential scaling of R 2 \ with the isotope neutron and proton num- 
bers was observed experimentally in the multifragmentation data from the 
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Figure 22: The isoscaling plots for the multifragmentation data from reac- 
tions of 124 Sn+ 124 Sn and 124 Xe+ 112 Sn at 28 A MeV [HI El at three different 
angular ranges. 



reactions of high-energy light particle with massive target nucleus |53| W2\ 
or from the reactions of mass symmetric projectile and target at intermedi- 
ate energies jM] and such behavior is called isoscaling [51] ( the parameters 
a,P,a',P' being called isoscaling parameters ). Isoscaling behavior was fur- 
ther reported in the heavy residue data [HI] and also in fission data jHH]- It 
was shown that the values of isoscaling parameters can be related to symme- 
try energy [2UE2], to the level of isospin equilibration [53] and to the values 
of transport coefficients [S3]. An example of isoscaling behavior of the frag- 
ments with Z=3-8 detected in the reactions of 124 Sn+ 124 Sn and 124 Xe+ 112 Sn 
at 28 AMeV [SHI EH] is given in Fig. [221 The isoscaling behavior is rather reg- 
ular at all three angular ranges, the isoscaling parameter a increases with the 
detection angle. Such a trend of the isoscaling parameter a can be possibly 
related to the fact that with increasing angle the fragments originate from 
still more damped collisions, with increasing level of isospin equilibration 
between projectile- and target-like fragment |56| . 

In Fig. [213 are presented the isoscaling data from statistical decay of 
hot quasiprojectiles from the reactions 28 Si+ 124,112 Sn at projectile energy 50 
AMeV [57J. Observed charged particles with Z<5 were isotopically resolved 
and total observed charge was close to the charge of the projectile ( Z=12- 
15 ) jI)El EH- Due to the restricted set of identified fragment species, the 
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Figure 23: The isoscaling plots from the reactions of 28 Si+ 124 ' 112 Sn at 50 
AMeV for the full set of isotopically resolved quasi-projectiles and for five 
bins of quasi-projectile excitation energy per nucleon |oT] . 



isoscaling plots are presented, according to formula (JT2J, as a function of 
isospin asymmetry. The isoscaling plots are presented not only for the full 
data but also for five bins of excitation energy. The slopes depend on the 
excitation energy in a similar way as the slope of the dependence of isobaric 
ratio Y( 3 H)/Y( 3 He) on the quasi-projectile N/Z observed in |Hj ( see Fig. EDI 
). When compared to similar results at 30 AMeV, the isoscaling parameters 
f3' for individual excitation energy bins do not depend on the projectile energy 
and the difference of the slopes of the inclusive data is caused by different 
excitation energy distributions of the hot quasi-projectiles at two projectile 
energies ( see Fig. 1131). 

As demonstrated, isoscaling appears to be a global feature of the multi- 
fragmentation data and the values of isoscaling parameters show sensitivity 
to properties of the hot source. It is of interest to explore such a sensitivity 
in the studies of thermodynamical properties of the multifragmenting source. 
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Thermo dynamical properties of the hot quasi- 
projectiles 

In the previous section we presented an overview of the methods by which 
the isotopic trends of the fragment data can be used to extract thermody- 
namical observables of the hot multifragmenting source. Here we present a 
systematic investigation of thermodynamical properties of the well character- 
ized hot quasi-projectiles with mass A=20-30, originating from the reactions 
28 Si+ 124 ' 112 Sn at projectile energy at 30 and 50 AMeV [TU]. Signatures and 
argumentation supporting the conclusion that multifragmentation of such 
quasiprojectiles is a statistical process will be presented. Furthermore, the 
signals for the isospin-asymmetric liquid-gas phase transition in the nuclear 
matter as a process underlying to multifragmentation will be investigated 
and implications for the nuclear equation of state at sub-saturation densities 
will be discussed. 

Production mechanism and equilibrium 

In order to establish whether the fragment data from reactions 28 Si+ 124,112 Sn 
at projectile energy 50 and 30 AMeV originate from statistical decay of 
hot quasiprojectiles, it is necessary to investigate the dynamical properties 
of the projectile-like source and to determine the production mechanism. 
The observed fragment data |581 ITU] provide full information ( with the ex- 
ception of emitted neutrons ) on the decay of thermally equilibrated hot 
quasi-projectiles with known mass ( A=20-30 ), charge, velocity and excita- 
tion energy. Detailed investigation of the reaction mechanism ^U] allowed to 
establish a dominant reaction scenario. Excellent description of the fragment 
observables was obtained using the model of deep-inelastic transfer ( DIT ) 
[ST?] for the early stage of collisions and the statistical multifragmentation 
model ( SMM ) [S] for de-excitation. The model describes well the dynam- 
ical properties of the reconstructed quasi-projectile such as center of mass 
velocity, excitation energy ( see Fig. ) and isospin-asymmetry ( see Fig. 
I2H). The fragment observables such as multiplicity, charge distributions and 
mean values of N/Z for a given charge were also reproduced reasonably well 
[TUj . Thus the data can be considered as well understood in terms of reaction 
mechanism. 

The observation that a wide range of observables is in agreement with the 
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Figure 24: Experimental ( solid circles ) and simulated ( solid lines ) mass 
distributions for the fully isotopically resolved quasiprojectiles with Z tot = 
14 HD|, (a) - 28 Si(30AMeV) + 112 Sn, (b) - 28 Si(30AMeV) + 124 Sn, (c) - 
28 Si(50AMcV) + 112 Sn : (d) - 28 Si(50AMcV) + 124 Sn. 

dominant reaction scenario of the deep inelastic transfer ( binary collision 
with intense nucleon exchange between the projectile and target nuclei, each 
of them represented by a statistically equilibrated Fermi gas ) followed by 
statistical multifragmentation ( assuming statistically equilibrated freeze-out 
configuration ) allows to conclude that the assumption of statistical equilib- 
rium is applicable and that thermodynamical observables can be defined for 
the hot quasi-projectile. The experimentally observed exponential scaling 
behavior shown in Figs. 1201 an d ESI further demonstrates compliance with 
grand-canonical picture. The contribution from non-equilibrium processes 
such as pre-equilibrium emission was shown to be weak pUj- The number 
of emitted neutrons, which are not detected, is, according to simulation, be- 
tween one and two per event and thus no significant distortion of the picture 
can be expected. The observed quasi-projectile multifragmentation data can 
thus be used for investigation of thermodynamical properties of the multi- 
fragmentation source and of the possible phase transition. 
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Thermodynamical observables and isospin-asymmetric 
phase transition 

In order to investigate thermodynamical properties of the hot multifragmen- 
tation source, one needs to characterize a set of observables, including exten- 
sive ones such as thermal excitation energy ( heat ) and intensive ones such 
as temperature and chemical potentials. 
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Figure 25: Dependence of the temperature T on excitation energy e* app . Solid 
squares - temperature determined using the double isotope ratio thermometer 
d,t/ 3 He, 4 He. Open squares - temperature determined from the dependence 
of the yield ratio Y( 3 H)/Y( 3 He) on the N/Z ratio of the quasiprojectile 

For the quasi-projectile data from the reactions 28 Si+ 124 ' 112 Sn at projec- 
tile energy 30 and 50 AMeV, temperature was determined both using dou- 
ble isotope ratio method and the N/Z- dependence of the isobaric yield ratio 
Y( 3 H)/Y( 3 He) 03]. The results are shown in Fig. |2ID Both methods seem to 
provide consistent results. Double isotope ratio temperature exhibits a short 
plateau starting around e* pp =4.5 AMeV, while the Y( 3 H)/Y( 3 He) tempera- 
ture seems to exhibit similar feature at somewhat higher excitation energy. 
The presence of a plateau in the caloric curve was interpreted in the litera- 
ture [nni as a signature of the first-order phase transition, nevertheless the 
influence of additional effects such as non-thermal flow or variation of the 
source size may affect the shape of the observed caloric curve |H3] and such 
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a signal can not be considered as an unambiguous signature without provid- 
ing further evidence. For the data in Fig. ESI the influence of non-thermal 
flow and variation of the source size can be excluded due to well understood 
reaction mechanism. 
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Figure 26: (a) The dependence of observable f3'T ( = A(jU n — ) on relative 
excitation energy of the quasiprojectile [HZj ( full squares ). Horizontal line 
shows estimate of (3'T obtained when assuming \i v ^ n w —S p , n . Dashed curve 
represents the expected trend for homogeneous system, (b) Corresponding 
caloric curve j33]. Vertical line indicates the position of turning point. 

Of the two methods for temperature determination, the Y( 3 H) /Y( 3 He) 
temperature relies on additional assumption concerning chemical potentials. 
The insight into behavior of the chemical potentials can be obtained using 
the isoscaling analysis. The values of the isoscaling parameter (3' for the 
data in Fig. ESI can be related, according to formula (fT2"|) . to the isovector 
component of the chemical potential. In Fig. EEH we present an estimate 
of the difference of isovector chemical potential using the observable f3'T, 
canceling out a trivial 1 /T-dependence of the isoscaling parameter . The 
double isotope ratio temperature [JS] from the isotopic ratios Y( 2 H)/Y( 3 H) 
and Y( 3 He)/Y( 4 He), which is independent of assumptions concerning chem- 
ical potentials, was used. The horizontal line represents an estimate of the 
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zero temperature value of A(/i n — by the proton and neutron separa- 
tion energies ( /i P) „ ~ — S P}Tl ) of reconstructed quasi-projectiles, using the 
known N/Z equilibration and the correction for neutron emission from back- 
tracing of DIT/SMM simulations [TO] - The zero temperature estimates of 
A(/i n — fMp) for subsets of data in all excitation energy bins are consistent 
with the horizontal line. Any experimental deviations from the horizontal 
line can be understood as a non-trivial dependence representing the details of 
de-excitation. The experimental dependence of f3'T indicates initial decrease 
at low excitation energies ( consistent with expansion of the homogeneous hot 
source ), further it exhibits a turning-point at 4 AMeV followed by increase 
toward the zero temperature value at 6-7 AMeV. The estimate of (3'T for 
homogeneous system ( dashed curve ) was obtained assuming the Fermi-gas 
p 2 / 3 -dependence and using the estimate of free volume by formula © for 
each excitation energy bin. The trend explains a reasonable success of the 
[J> P ,n ~ — S p , n approximation used to extract the mirror nucleus temperature 
in Fig. |2I3 where the viability of the approximation at high temperatures 
was supported by a possible counterbalance of two effects, namely expansion 
of the hot nucleus ( leading to a decrease of the absolute values of chemical 
potentials ) and chemical separation into an isospin symmetric heavy fraction 
and an isospin asymmetric nucleon gas, leading to increase of \i n — \i v of the 
dilute phase ( nucleon gas ). Since the isoscaling parameter depends directly 
on the free nucleon chemical potentials, the turning-point at 4 AMeV can be 
understood cLS cL si gnal of the onset of chemical separation which reverts the 
decrease of the free nucleon chemical potential consistent with expansion of 
the homogeneous system. As demonstrated by comparison with the caloric 
curve in Fig. EoT o ( obtained using double isotope ratio thermometer without 
assumption concerning chemical potentials ) the onset of chemical separation 
is correlated to the onset of the plateau in the caloric curve, thus signaling 
that chemical separation is accompanied by latent heat. Such a behavior can 
be expected in the first-order phase transition. The dotted lines in Fig. l2B*b 
indicate the discontinuity in the heat capacity of the system. 

The behavior, observed in Fig. [20] is consistent with the first-order phase 
transition. Several signals of the first-order phase transition have been re- 
ported in the recent experimental studies, such as negative micro-canonical 
heat capacity |62| . characterized by abnormal kinetic energy fluctuation, and 
the "fossil" signal of spinodal decomposition [02], characterized by abnormal 
production of equally sized fragments. These indications of the first-order 
phase transition are based on careful event-by-event analysis of the multifrag- 
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Figure 27: (a) - Multiplicity of charged fragments (squares), LCPs (cir- 
cles) and IMFs (triangles) versus N/Zgp. Corresponding lines represent 
the DIT/SMM calculation, (b) - Experimental ratio of the mean values 
of N/Z of LCPs and IMFs versus the N/Zgp. Data are given for the reaction 
28 Si(50AMeV)+ 112 Sn 



mentation data combined with model assumptions and results of theoretical 
simulations. The isospin-asymmetry of the fragments is not explored in such 
analyzes. The behavior in Fig. |2Sl relates the signal of chemical separa- 
tion to the plateau of the caloric curve, thus providing a signal of isospin- 
asymmetric first-order phase transition. Further support for such conclusion 
can be deduced from Fig. l2Tb . demonstrating the inhomogeneous distribu- 
tion of isospin-asymmetry among dense and dilute phases represented by 
light charged particles and intermediate mass fragments, respectively. The 
effect is strongly dependent on the N/Z ratio of the quasiprojectile. So are 
multiplicities of LCPs and IMFs in Fig. l2Tk which further indicate that with 
increasing isospin-asymmetry of the quasiprojectile the dilute phase becomes 
more isosp in- asymmetric and abundant. The observations related to Figs. 
l2*o1 1771 are consistent with theoretical predictions of the isospin-asymmetric 
phase transition in nuclear media. The estimated density of the homogeneous 
system at the turning point in Fig. l2*Uk is consistent with theoretical esti- 
mates of the position and shape of the spinodal contour, typically at almost 
constant total density 0.6po in a wide range of asymmetries jjj. Thus, in- 
side the spinodal region the homogeneous nuclear medium is quickly replaced 
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by both isospin and spatially inhomogeneous system. Such a conclusion is 
further supported by a comparison of the temperature obtained using the 
double isotope ratio thermometer d,t/ 3 He, 4 He jHj with kinematic tempera- 
tures of p,d,t, 3 He, 4 He obtained from Maxwellian fits of the particle kinetic 
energy spectra in the quasiprojectile frame ( see Fig. |2HI )• While at low 
excitation energies the double isotope ratio thermometer d,t/ 3 He, 4 He rep- 
resents essentially an average of the kinematic temperatures of individual 
species, in the plateau region occurs a transition to the regime where the 
double isotope ratio temperature represents the kinematic temperature of 
protons ( and possibly of the nucleon gas ). Such a transition, just above the 
turning-point where the isovector part of chemical potential deviates from 
the trend of homogeneous system (as shown in Fig. suggests that the 
isotopic composition of fragments is in equilibrium with the free nucleon gas 
what is a basic assumption used in grand-canonical models such as SMM |5 S . 
The higher values of kinematic temperatures of d,t, 3 He, 4 He suggest that the 
clusters are created prior to such equilibration ( possibly by spinodal decom- 
position ). On the other hand, the low energy behavior can be explained 
by an assumption that the isotopic composition is determined at the instant 
when clusters are emitted what is in good agreement with our understanding 
of de-excitation of the compound nucleus [T2] . The two different values of 
heat capacities in low- and high-energy parts of the caloric curve, suggested 
in Fig. |2Sb, can thus be understood as representing a homogeneous nuclear 
liquid and an isospin-asymmetric free nucleon gas, respectively. 

Apart from signatures of the first-order phase transition, several experi- 
mental works [Hll EI, exploring observed mass ( charge ) distributions, report 
a critical behavior in the multifragmentation data, which can be understood 
as a signal of the second-order phase transition at the critical temperatures 
between 5-10 MeV [US]- However, multifragmentation studies employing 
the statistical model analysis of charge distribution suggest the value 
of critical temperature above 15 MeV. Besides of raising obvious questions 
concerning the applicability of used methods of analysis, such contradictory 
observations may in principle suggest a complex phase diagram, possibly 
analogous to the Ising model, where both first-order and second-order can 
be observed. It is also possible that the observed critical behavior is a spe- 
cific feature of the finite system, as suggested in theory \§7\. The recent 
experimental work where the critical behavior is related to the size of largest 
fragment as an order parameter jHHl supports such assumptions. The finite- 
size effect can play role also in the multifragmentation of quasiprojectiles 
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Figure 28: Comparison of temperature obtained using the double iso- 
tope ratio thermometer d,t/ 3 He, 4 He |44] with kinematic temperatures of 
p,d,t, 3 He, 4 He in the cm. frame of the quasi-projectiles. 

with A=20-30 and thus similar analysis for heavier sources is necessary to 
evaluate the role of such effect and establish the link to the phase transition 
in the isospin-asymmetric nuclear matter at sub-saturation densities. It may 
be however noted that theoretical studies of small systems suggest that 
van der Waals-like loops observed in infinite systems are affected in transi- 
tion to finite system only by small terms and thus the conclusions for finite 
systems can be in principle applicable also for infinite case. 

Summary and conclusions 

In this article we presented an overview of the recent progress in the studies 
of nuclear multifragmentation. Theoretical and experimental concepts rele- 
vant to investigation of isotopic trends in nuclear multifragmentation were 
described and the possibilities to extract the physical information on the nu- 
clear equation of state were outlined. A special emphasis was put on the 
methods how the isotopic composition of fragments can be used to extract 
the values of thermodynamical observables of the system undergoing mul- 
tifragmentation such as temperature and chemical potentials ( free nucleon 
densities ). Various methods for extraction of thermodynamical variables 
were presented and possibilities of their use were documented on examples 
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of multifragmentation data. It was shown that the yield ratios of specific 
isotopes from statistical multifragmentation typically exhibit an exponential 
scaling with various observables. The logarithmic slopes of such dependences 
can be related to thermodynamical observables such as grand-canonical tem- 
peratures ( isotopic thermometers ) or chemical potentials ( isoscaling ). 
Finally, an extensive overview of investigations of multifragmentation of the 
fully reconstructed quasi-projectiles is presented. The dominant reaction 
mechanism is determined as binary nucleus-nucleus collisions with intense 
dissipation and heating via nucleon exchange. The implications of the frag- 
ment data concerning the isospin-asymmetric liquid-gas phase transition in 
the nuclear matter are discussed. Multiple correlated signals point toward 
observation of a first-order phase-transition in the quasi-projectile multifrag- 
mentation data. 
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